Topologically trapped vortex molecules in Bose-Einstein condensates 
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In a numerical experiment based on Gross-Pitaevskii formalism, we demonstrate unique topolog- 
ical quantum coherence in optically trapped Bose-Einstein condensates (BECs). Exploring the fact 
that vortices in rotating BEG can be pinned by a geometric arrangement of laser beams, we show 
the parameter range in which vortex- antivortex molecules or multiquantum vortices are formed as a 
consequence of the optically imposed symmetry. Being low-energy states, we discuss the conditions 
for spontaneous nucleation of these unique molecules and their direct experimental observation, and 
provoke the potential use of the phase print of an antivortex or a multiquantum vortex when realized 
in unconventional circumstances. 

PACS numbers: 03.75.Lm, 67.85. Hj, 47.32. cd 



During the last decade, intense experimental and the- 
oretical activity has induced dramatic achievements in 
Bose-Einstein condensation in trapped alkaU-metal gases 
at ultralow temperatures [H. The atomic Bose-Einstein 
condensates (BECs) differ fundamentally from the hc- 
Uum BEC, having significantly non-uniform density, very 
high low-temperature condensate fraction, and being 
pure and dilute, so that interactions can be accurately 
parametrized in terms of the scattering length. As a re- 
sult, a relatively simple nonlinear Schrodinger equation 
(the Gross-Pitaevskii (GP) equation) gives a precise de- 
scription of the atomic condensates and their dynamics 

i- 

It is already well established that the GP equation has 
rotating solutions with a well-defined angular momen- 
tum, defining the vortex states Vortices as topolog- 
ical singularities in BEC have been readily observed in 
experiment f^, and are also found in liquid helium Q 
and in superconductors Q. The exciting new property 
of atomic BECs is that intrinsic interactions in the sys- 
tem are directly dependent on the mean particle density 
01 , which gives them a unique degree of freedom for stud- 
ies of vortex matter. Another interesting aspect of dilute 
Bose condensates is that their density can be locally sup- 
pressed optically. As experimentally shown by Tung et 
al. [1], vortex lines in BEC show great affinity towards 
weak spots created by focussed laser beams, which re- 
sults in a very profound vortex pinning and provides yet 
another tool for vortex manipulation [9|. 

Vortices in BEC are accompanied by the appropriate 
distribution of the phase ip of the order parameter ijj, 
where the number of successive 27r changes along the 
perimeter corresponds to the vorticity of the vortex state, 
i.e. the topological charge. Therefore, using the known 
phase distribution, particular vortex states can be prede- 
fined in the system by the so-called phase imprint tech- 
nique [lo| . where an uniform light pulse projects a de- 
signed mask onto the condensate. However, the gradient 
of the imprinted phase (S/ip) must always be parallel to 



the rotational velocity of the condensate, i.e. the an- 
tivortex is fundamentally unstable for a positive angu- 
lar momentum. Note also that imprinted multiple 2tt 
change, i.e. a multiquantum vortex, was never found sta- 
ble against decomposition into single vortices in conven- 
tional BEC setups [nj. 

Motivated by these recent developments in vortex 
physics, in this paper we exploit a uniquely defined BEC 
system that combines above essentials; our atomic BEC 
is confined, exposed to engineered spatial pinning, and is 
cylindrically rotated. Knowing that vortex solitons are 
strongly influenced by symmetry in non-linear media , 
we impose a particular symmetry on the system by cho- 
sen arrangement of the vortex-trapping optical beams. 
As we will show, this results in the nucleation of vor- 
tex molecules of matching symmetry in the ground state, 
which may comprise unconventional phase profiles con- 
taining a stable antivortex, or a multiquantum vortex. 

For the purpose of this work, we solve the stationary 
Gross-Pitaevskii (GP) equation, in a rotating frame with 
frequency around the z-axis: 



- — A + V, + g\lP\'-nL,,-^l]i,^0, (1) 



where Vc stands for confinement potential, g = ^■nh^a/m 
is the non-linearity with a being the s-wave scattering 
length, Lz = —ih(xdy — ydx) is the angular momentum 
operator, /i denotes chemical potential, and ip is the or- 
der parameter normalized to the number of particles in 
the condensate N . The energy of a particular state is 
then given hy E = ji — ^ j We take the usual 

parabolic-like confinement potential characterized by the 
frequency uj±^ in the {x, y) plane and lo^ in the z-direction. 
Our numerical method treats the GP equation in all three 
dimensions, but in the present work we take lOz ^ t^_L 
forcing the condensate in a quasi two-dimensional, oblate 
(pancake) shape. In that case, Eq. ([T]) can be written in 
dimensionless and discretized form as 
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where 6; 
grid 



and the 



are the lattice constants of the Cartesian 



link variable is defined as 



r 7"m,n 
^ OL—x^y 
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with A,^(y) = {—)^VLy{x). The 
distances, angular velocity and energy are expressed in 
the fundamental scales of the harmonic trap, i.e. = 
and Ef) ~ fkoj 

2aN. 



2D nonlinearity g' 



respectively. The 



l2iTmuiz/h (found by averag- 
ing the GP equation over the z-direction) now depends 
on N since we chose wave function normalized to unity. 
Therefore, for fixed 5, is a measure of the importance 
of interactions in the condensate, and directly reflects 
on the vortex phase diagram. Note also that for low 
non-linearities, i.e. when g' is small, vortex states in 
parabolic confinement are stable only for velocities of 
the condensate close to f2 = wj,, i.e. when the gas dis- 
solves due to centrifugal forces. To enhance the vortex 
stability, we have chosen a modified confinement poten- 
tial Vc = -|- j^r^. Such empowered confinement en- 
ables the condensate to survive higher angular velocities 
(f2 > OJl_)^ and is feasible in experiment by combining 
optical and magnetic trapping (e.g. for ^^Rb, a magnetic 
trap with frequency 75.5 x 27r Hz in combination with a 
laser beam with waist 25 ^m and power 1.2 mW [l3|). 

Vortex- antivortex molecules as the ground state. In 
this paper, we consider '^^Na. BEC, with somewhat lower 
interaction than ^^Rb, and consequently lower non- 
linearity. With the above expression for Vc, we use 
the effective perpendicular parabolic confinement with 
uj± = 4.39 X 2tt Hz, and the unit of length becomes 
ro = 10 /im. To obtain the pancake shape of the conden- 
sate, we take uj^ = lOwj^. The non-linearity parameter 
for a condensate consisting of 10^ sodium atoms becomes 
g' = 44. In what follows, we implement laser pinning in 
the formalism, and model it as a Gaussian potential well 
Vp = exp(— 16(r/wp)2) (r is the distance from the pinning 
center) as depicted schematically in the inset of Fig. [Tfc. 
We use a set of Np laser beams on a single ring, and 
compare first the energy in the absence and presence of 
pinning, as shown in Fig. [TJ:. The resulting total po- 
tential in the BEC, introduced by confinement and laser 
beams, is taken as 



V = Vcir^ + 



with -fj = dp[cos{j2TT /Np)ex + s\n[j2T: /Np)ey]. With 
increasing angular velocity lo = il/tjj^, vortices stabi- 
lize in the system, and each vortex state exhibits a dis- 
tinct energy. In the Np = case, for chosen parame- 



ters, individual vortices form multi-vortex patterns de- 
noted by {Lc,L), where Lc gives the number of vor- 
tices in the central part of the condensate i.e. Lc = 
lim f d (log 1^ j , and L is the total vorticity. In 

the Np ^ case, the ordering of laser beams imposes its 
symmetry on the vortex states, with Np vortices pinned 
by the laser beams, and the remaining ones sitting in the 
central region of the condensate. As a remarkable phe- 
nomenon, the system may preserve the symmetry even 
for L < Np-, e.g. for L = A^p — 1, still Np vortices nu- 
cleate and are pinned, but must be accompanied by a 
central antivortex in order to simultaneously match the 
total angular momentum and the symmetry of the pin- 
ning potential. Energy levels of such vortex-antivortex 
molecules are indicated in Fig. [1] for Np = 3 and Np ~ 5 
by the dotted (blue) curves. 

Vortex-antivortex (V-Av) phenomena has already been 
of interest in BEC, mainly due to the fact that the exci- 
tation energy of a vortex-antivortex pair in non-rotating 
condensate is lower than the one of a single vortex. For 
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FIG. 1: (Color online) Energy versus angular velocity of a 
^^Na condensate with g'=20 (i.e. with A'' « 5000 atoms), in 
the case of no laser pinning (A^'p = 0, a), and three {Np = 3, 
b), and five {Np = 5, c) laser beams in a polygonal setting 
(with dp = 0.5, Wp — 0.4, Up = 500). Different lines cor- 
respond to states with different total vorticity L. The colors 
and style of the lines indicate the vorticity in the center of the 
condensate (Lc). The inset depicts the used energy profile for 
each pinning site, and notation of the parameters. 
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that reason, it was assumed that the actual nucleation 
of a vortex occurs through the generation of a V-Av 
pair and the subsequent expulsion of the antivortex 14 1 . 
While this may hold in ideally non- interacting BEC, V- 
Av appearance is seriously hampered by even small per- 
turbations in the system, particularly in the rotating one. 
Yet, Crasovan et al. demonstrated that a V-Av pair may 
be stable in an interacting non-rotating BEC, however 
only as an excited state (ground state remains vortex- 
free) [3l and the dynamics of such a state was recently 
analyzed in Ref. [16[ . In addition, vortex-antivortex clus- 
ters were studied in Ref. , the stabilization of vortex- 
antivortex lattices was considered in Ref. [3l and a su- 
perposed vortex-antivortex state was investigated in Ref. 
1^ . however all only for the case of non-rotating systeuia. 
Here we show that states comprising an antivortex can 
be stabilized as the ground state in a rotating BEC by 
engineered pinning. Moreover, by using a different num- 
ber of laser beams, antivortex can be realized at different 
angular velocities. Fig. [2] shows the w-stability range of 
V-Av states with one antivortex found for Np = 3 — 7, 
thus all states having vorticity L = Np — 1. With increas- 
ing Np, not only V-Av states appear for higher angular 
momenta, but they can also be spread over a larger area 
of the condensate. When artificially spreading the laser 
beams further apart (increasing dp, see Fig. [1]), vortices 
spontaneously follow and increase their distance from the 
central antivortex (dyav ) which facilitates the experimen- 
tal verification of the V-Av molecule. Fig. [2] also shows 
the maximal dyav that can be reached by increasing dp 



for 



0.5, Op = 500 and g' = 40. The insets of Fig. [2] 




FIG. 2: (Color online) Histogram of vortex-antivortex (V- 
Av) distance in stable molecules found for different number 
of equidistant laser beams on a ring, as a function of angular 
velocity for g' = 40. Maximal V-Av separation increases for 
more used pinning beams, and insets show two oblique views 
of the superimposed density/phase 3D plots of the condensate 
for u> = 1.71 and iV„ = 7. 



show isosurfaces of density 10 ^ on which the phase dis- 
tribution is superimposed for a condensate in the V-Av 
state for Np = 7, g' = 22 and lj = 1.71. 

Novel phase transitions. As much as V-Av states are 
novel and exciting study objects, our system actually ex- 
hibits even richer vortex structures beyond V-Av phe- 
nomena. To illustrate this, we constructed a full vortex 
phase diagram, shown in Fig. [3] as a function of uj and 
non-linearity g' , for Np = 4, i.e. a square arrangement of 
laser beams. As shown before, for constant g' , the vor- 
ticity increases with increasing lo. Note however that our 
strategic pinning setup enforces the A^p-fold symmetry of 
the vortex states, which leads to specific transitions in 
the central region of the condensate. As illustrated by 
a color gradient in Fig. [3l the central vorticity changes 
as Lc = mod (L + I, Np) — 1, where an antivortex may 
nucleate as discussed above. Nevertheless, iVp-fold sym- 
metry tends to remain preserved for other vorticities as 
well, which results in the compression of excess L — kNp 
[k G N'^) vortices in the center of the condensate into a 
multiquantum vortex (also called giant vortex). This re- 
sult complements earlier studies of such vortices in con- 
ventional BEC setups, where a multi-quanta vortex was 
found unstable towards decay into single-quantum vor- 
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FIG. 3: (Color online) Equilibrium vortex phase diagram as 
a function of the angular velocity and non-linearity of the 
condensate, for Np = 4, dp = 0.71, Wp = 0.4 and Up = 500. 
Coloring scheme illuminates the vorticity in the center of the 
condensate, between the pinning beams. Solid lines show the 
transitions in the ground state, while dashed curves illustrate 
the whole region of stability for different vortex states. In 
white areas, the vortex structure does not fully obey the four- 
fold symmetry of the pinning. 
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tices For Np = 4, we show the 0,1,2,-1,0,1,2, -1 
sequence with increasing lu in Fig. [3l for total vorticity 
L = — 7. Doubly quantized vortex was realized for 
L = 2 and L = 6 states in the center of the condensate. 

While the existence of stable multiquantum vortices 
was already predicted in several theoretical papers for 
condensates in a confinement potential with a quartic 
component (20| . our giant vortices are solely a conse- 
quence of the symmetry induced by the pinning centers, 
as can be deduced from Fig. \M^d) for the L = 6 = 4 + 2 
state: The giant vortex state is not stable when Wp is too 
small. By doing two separate simulations of a conden- 
sate in conventional parabolic confinement - thus with- 
out quartic component - exposed to 4 pinning beams, we 
were able to confirm this statement: We found a giant 
vortex in both the ground states for L = 2 (with g = 20, 
Lu = 0.85, dp = 0.8, Wp = 0.6, ap = 50) and L = 6 = 4 + 2 
(g = 40, w = 0.85, dp = 1.3, Wp = 1.6, Up = 200). This 
general behavior holds for arbitrary number of pinning 
beams. In other words, vortex molecules comprising mul- 
tiquantum vortices are stable in our system, contrary to 
conventional behavior of BEC. Moreover, we report here 
a unique way of engineering locally the phase imprint of 
the condensate, by either changing angular velocity, or 
the non-linearity (e.g. by condensing more constituent 
atoms). The limits of this procedure are: (i) the decon- 
finement of the condensate at large lu, and (ii) the maxi- 
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FIG. 4: (Color online) Influence of the optical parameters of 
the laser beams on the stability of the special vortex molecules 
found for Np — 4. (a-c) V-Av stability and spacing in a two- 
parameter space, thus for either fixed amplitude, separation, 
or focused width of the beams {g' = 40, u; — 1.25). (d) The 
stability region of the doubly quantized vortex in the center 
of a L = 6 state, for different values of ap {g' = 40, oj = 1.6). 



mal g' at which vortex- vortex interaction overwhelms the 
optical pinning, resulting in a vortex structure that does 
not fully obey the imposed geometry (white areas in Fig. 

ED. 

Stability and control of special vortex states. Obvi- 
ously, the exact parameters of the strategically placed 
laser beams are very deterministic for the resulting vortex 
structure in the condensate. Here we briefly discuss the 
influence of beam intensity, distance between adjacent 
beams, and their focussed waist on V-Av and multiquan- 
tum states. Our results are summoned in Fig. |4l Figs. 
(a)-(c) relate to V-Av molecules in fourfold pinning, and 
help visualize their stability and size in a two-parameter 
space. As a general conclusion, we find that a larger am- 
plitude of the pinning (ap, see inset of Fig. [T]) is always 
beneficial to the V-Av state. Simply, stronger imposed 
pinning reflects better its symmetry on the vortex state. 
For the same reason, strong overlap between pinning pro- 
files of the beams should be avoided (i.e. Wp < 3dp is 
desirable) . As shown in Fig. 01 the distance between the 
central antivortex and pinning vortices can be increased 
by increasing the spacing between the beams (dp); this is 
nevertheless a limited option, as vortices may not follow 
the imposed pinning at large distances from the center of 
the condensate. For considered parameters, we realized 
maximal V-Av distance of approximately O.Qrg, which is 
of the order of lO^m for the sodium condensate, and thus 
feasible for experimental observation. The used beam 
width of Wp = 0.2 — 1.5ro and sub-rg distances for beam 
ordering are also very realistic values for the available 
experimental techniques. 

Fig. md) shows the stability region of a doubly quan- 
tized vortex in a fourfold geometry of the L = 6 state 
(4 X 1 + 2). Similarly to the V-Av case, imposed sym- 
metry of the pinning is essential for the nucleation of 
the double vortex. Relatively weaker pinning by wider 
laser beams can significantly enhance the stability of the 
multiquantum vortex, while the clarity of its separation 
from the surrounding vortices is best achieved for high 
beam amplitudes and narrow beam waists. However, the 
needed dp is quite low, as confinement between the pin- 
ning beams must remain firm for the stabilization of the 
multiquantum vortex. Such small spacing is unfavorable 
for experimental observation, but can be significantly im- 
proved in the case of larger Np and larger vorticity. This 
analysis and the complete investigation of all stable vor- 
tex states in this phenomenologically rich system will be 
presented elsewhere. 

Conclusions. In summary, using a novel concept of 
polygonal optical pinning, we demonstrated yet unpre- 
dicted vortex states in rotating atomic Bose-Einstein con- 
densates. We realized a stable antivortex for positive 
angular momentum, as well as a multiquantum vortex, 
which up to now was assumed to be unstable in BEC. 
Both are found as the ground state for a wide range of 
parameters. New ground-state transitions found by in- 
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creasing angular velocity of the condensate, by changes 
in the optical pinning setup, or by changing the number 
of constituent atoms, open up ways to further studies em- 
ploying quantum phase engineering and submicron 
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coherence and matter-wave interference effects 

The authors thank J. Tempere for a critical reading of 
the manuscript. This work was supported by the Flemish 
Science Foundation (FWO-Vl). M.V.M. acknowledges 
support from the EU Marie-Curie Intra-European pro- 
gram. The work was partially done during visit of F.M.P. 
to the Institute for Mathematical Sciences, National Uni- 
versity of Singapore in 2007. 



[5] 
[6] 
[7] 
[8] 
[9] 



Electronic address: |francois.peeters@ua.ac.be| 
M.H. Anderson, J.RrEnsher, M.H. Matthews, C.E. Wie- 
man, and E.A. Cornell, Science 269, 198 (1995); K.B. 
Davis, M.-O. Mewes, M.R. Andrews, N.J. van Druten, 

D. S. Durfee, D.M. Kurn, and W. Ketterle, Phys. Rev. 
Lett. 75 3969 (1995). 

E. M. Lifshitz and L.P. Pitaevskii Statistical Physics, Vol. 
2 (Pergamon, Oxford, 1980). 

G. Baym and C.J. Pethick, Phys. Rev. Lett. 76, 6 (1996). 
M.R. Matthews, B.P. Anderson, P.C. Haljan, D.S. Hall, 

C. E. Wieman, and E.A. Cornell, Phys. Rev. Lett. 83 
2498 (1999); M.W. Zwierlein, J.R. Abo-Shaeer, A. Schi- 
rotzek, C.H. Schunck, and W. Ketterle, Nature 435, 1047 
(2005). 

D. R. Tilley and J. Tilley, Superfluidity and Superconduc- 
tivity (Hilger, Bristol, 1986). 

M. Tinkham, Introduction to Superconductivity 
(McGraw-Hill, New York, 1996). 

F. Dalfovo, S. Giorgini, L.P. Pitaevskii, and S. Stringari, 
Rev. Mod. Phys. 71, 463 (1999) 

S. Tung, V. Schweikhard, and E.A. Cornell, Phys. Rev. 
Lett. 97, 240402 (2006). 

K. Kasamatsu and M. Tsubota, Phys. Rev. Lett. 97, 
240404 (2006). 



[10] J. Denschlag, J.E. Simsarian, D.L. Feder, C.W. Clark, 
L.A. Collins, J. CubizoUes, L. Deng, E.W. Hagley, K. 
Helmerson, W.P. Reinhardt, S.L. Rolston, B.I. Schneider, 
and W.D. Philips, Science 287, 97 (2000). 

[11] D.A. Butts and D.S. Rokhsar, Nature 397, 327 (1999); 
A.E. Leanhardt, A. Gorlitz, A. P. Chikkatur, D. Kielpin- 
ski, Y. Shin. D.E. Pritchard, and W. Ketterle, Phys. Rev. 
Lett. 89, 190403 (2002). 

[12] Y.V. Kartashov, A. Ferrando, A. A. Egorov, L. Torner, 
Phys. Rev. Lett. 95, 123902 (2005); V.M. Perez-Garcfa, 
M.A. Garci'a-March, and A. Ferrando, Phys. Rev. A 
75, 033618 (2007); M.A. Garcfa-M arch, A. Ferrand o, 
M. Zacares, J. Vijande, L.D. Carr, larXiv:0806.0947k l 
(2008). 

[13] V. Bretin, S. Stock, Y. Seurin, and J. Dalibard, Phys. 
Rev. Lett. 92, 050403 (2004). 

[14] C. Raman, J. R. Abo-Shaeer, J. M. Vogels, K. Xu, and 
W. Ketterle, Phys. Rev. Lett. 87 210402 (2001); Q. Zhou 
and H. Zhai, Phys. Rev. A 70, 043619 (2004). 

[15] L.-C. Crasovan, G. Molina- Terriza, J. P. Torres, L. 
Torner, V.M. Perez- Garcia, and D. Mihalache, Phys. 
Rev. E 66, 036612 (2002); L.-C. Crasovan, V. Veksler- 
chik, V.M. Perez-Garcfa, J. P. Torres, D. Mihalache, and 
L. Torner, Phys. Rev. A 68, 063609 (2003). 

[16] W. Li, M. Haque, and S. Komineas, Phys. Rev. A 77, 
053610 (2008). 

[17] M. Mottonen, S.M.M. Virtanen, T. Isoshima, and M.M. 

Salomaa, Phys. Rev. A 71, 033626 (2005); V. Pietila, M. 

Mottonen, T. Isoshima, J.A.M. Muhtamaki, and S.M.M. 

Virtanen, Phys. Rev. A 74, 023603 (2006). 
[18] G. Ruben, D.M. Paganin, and M.J. Morgan, Phys. Rev. 

A 78, 013631 (2008). 
[19] M. Liu, L.H. Wen, H.W. Xiong, and M.S. Zhan, Phys. 

Rev A 73, 063620 (2006) 
[20] A. Aftalion, and I. Danaila, Phys. Rev. A 69, 033608 

(2004); I. Danaila, Phys. Rev. A 72, 013605 (2005); A.L. 

Fetter, B. Jackson, and S. Stringari, Phys. Rev. A 71, 

013605 (2005); K. Kasamatsu and M. Tsubota, Phys. 

Rev. Lett. 97, 240404 (2006). 
[21] P.R. Berman, Ed., Atom Interferometry (Academic 

Press, Cambridge, 1997). 



